Improved third‐order weighted essentially nonoscillatory scheme

A new third‐order WENO scheme is proposed to achieve the desired order of convergence at the critical points for scalar hyperbolic equations. A new reference smoothness indicator is introduced, which satisfies the sufficient condition on the weights for the third‐order convergence. Following the truncation error analysis, we have shown that the proposed scheme achieves the desired order accurate for smooth solutions with arbitrary number of vanishing derivatives if the parameter ε satisfies certain conditions. We have made a comparative study of the proposed scheme with the existing schemes such as WENO‐JS, WENO‐Z, and WENO‐N3 through different numerical examples. The result shows that the proposed scheme (WENO‐MN3) achieves better performance than these schemes.     

Maximal Subrings and Covering Numbers of Finite Semisimple Rings

We classify the maximal subrings of the ring of n×n matrices over a finite field, and show that these subrings may be divided into three types. We also describe all of the maximal subrings of a finite semisimple ring, and categorize them into two classes. As an application of these results, we calculate the covering number of a finite semisimple ring.     

Domain decomposition methods in a geometrical multiscale domain using finite volume schemes

The heat equation is solved by using a finite volume discretization in a domain that consists of a two-dimensional central node and several one-dimensional outgoing branches. Several interface connection options to match the submodels set on the node and on the branches, with or without continuity, are looked at. For each of them, a monolithic scheme is defined, and existence and uniqueness of the solution is proved. New schemes are deduced, which are obtained through domain decomposition methods in the form of interface systems, with one or two unknowns per interface. A comparative systematic study is carried out from an algebraic and numerical point of view according to the interface conditions: Dirichlet, Neumann, or Robin. An efficient diagonal preconditioning is proposed.     

Time-dependent Ginzburg–Landau equations for multi-gap superconductors

We numerically solve the time-dependent Ginzburg–Landau equations for two-gap superconductors using the finite-element technique. The real-time simulation shows that at low magnetic field, the vortices in small-size samples tend to form clusters or other disorder structures. When the sample size is large, stripes appear in the pattern. These results are in good agreement with the previous experimental observations of the intriguing anomalous vortex pattern, providing a reliable theoretical basis for the future applications of multi-gap superconductors.     

Identifying the viscoelastic properties of soft matter from the indentation response of a hard film-soft substrate system

The indentation technique is widely used in measuring the mechanical properties of soft matter at the microscale or nanoscale,but still faces challenges by these unique properties as well as the consequent strong surface adhesion, including the strong nonlinear effect, unclear judgment of the contact point, difficulties in estimating the contact area, and the risk of the indenter piercing the sample. Here we propose a two-step method to solve these problems: lay a hard film on a soft matter, and obtain the viscoelastic properties of this soft matter through the indentation response of this composite structure. We first establish a theoretical indentation model of the hard film-soft substrate system based on the theory of plates, elastic-viscoelastic correspondence principle and Boltzmann superposition principle. To verify the correctness of this method, we measure the mechanical properties of the methyl vinyl silicone rubber(MVSR) covered by a Cu nanofilm. Finally, we test the effectiveness and error sensitivity of this method with the finite element method(FEM). The results show that our method can accurately measure the mechanical properties of soft matter, while effectively circumventing the problems of the traditional indentation technique.     

Finite element computations of viscoelastic two-phase flows using local projection stabilization

A three-field local projection stabilized (LPS) finite element method is developed for computations of a three-dimensional axisymmetric buoyancy driven liquid drop rising in a liquid column where one of the liquid is viscoelastic. The two-phase flow is described by the time-dependent incompressible Navier-Stokes equations, whereas the viscoelasticity is modeled by the Giesekus constitutive equation in a time-dependent domain. The arbitrary Lagrangian-Eulerian (ALE) formulation with finite elements is used to solve the governing equations in the time-dependent domain. Interface-resolved moving meshes in ALE allows to incorporate the interfacial tension force and jumps in the material parameters accurately. A one-level LPS based on an enriched approximation space and a discontinuous projection space is used to stabilize the numerical scheme. A comprehensive numerical investigation is performed for a Newtonian drop rising in a viscoelastic fluid column and a viscoelastic drop rising in a Newtonian fluid column. The influence of the viscosity ratio, Newtonian solvent ratio, Giesekus mobility factor, and the Eötvös number on the drop dynamics are analyzed. The numerical study shows that beyond a critical Capillary number, a Newtonian drop rising in a viscoelastic fluid column experiences an extended trailing edge with a cusp-like shape and also exhibits a negative wake phenomena. However, a viscoelastic drop rising in a Newtonian fluid column develops an indentation around the rear stagnation point with a dimpled shape.     

A framework to perform high-order deconvolution for finite-volume method on simplicial meshes

In this article, a new framework to design high-order approximations in the context of node-centered finite volumes on simplicial meshes is proposed. The major novelty of this method is that it relies on very simple and compact differential operators, which is a critical point to achieve good performances in the High-performance computing context. This method is based on deconvolution between nodal and volume-average values, which can be conducted to any order. The interest of the new method is illustrated through three different applications: mesh-to-mesh interpolation, levelset curvature computation, and numerical scheme for convection. Higher order can also be achieved within the present framework by introducing high-rank tensors. Although these tensors feature much symmetries, their manipulation can quickly become an overwhelming task. For this reason and without loss of generality, the present articles are limited to third-order expansion. This method, although tightly connected to the k-exact schemes theory, does not rely on successive corrections: the high-order property is obtained in a single operation, which makes them more attractive in terms of performances.     

L1-type smoothness indicators based weighted essentially nonoscillatory scheme for Hamilton-Jacobi equations

This article presents an improved fifth-order finite difference weighted essentially nonoscillatory (WENO) scheme to solve Hamilton-Jacobi equations. A new type of nonlinear weights is introduced with the construction of local smoothness indicators on each local stencil that are measured with the help of generalized undivided differences in L¹-norm. A novel global smoothness measurement is also constructed with the help of local measurements from its linear combination. Numerical experiments are conducted in one- and two-dimensions to demonstrate the performance enhancement, resolution power, numerical accuracy for the proposed scheme, and compared it with the classical WENO scheme.     

基于内聚力模型的快速荷载下CFRP-混凝土界面剪切剥离模拟

鉴于纤维增强复合材料(Fiber Reinforced Polymer, FRP)已在混凝土结构加固中普遍应用, 能有效提高结构的承载和变形能力. 为揭示动载作用下FRP-混凝土界面的剥离机理, 本文基于ABAQUS平台, 采用内聚力模型模拟CFRP-混凝土界面层, 实现了快速荷载下CFRP-混凝土界面剥离的高效模拟. 结果表明: CFRP表面应变在加载过程中由加载端向自由端方向传递, 随着加载速率的提高, 界面承载能力也随之提高; 界面峰值剪应力也存在显著的应变率效应; 模拟结果与试验结果基本吻合, 说明了新方法的有效性.